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SUMMARY 

Droplet growth equations are reviewed in the free-molecular, transi- 
tion, and continuum flow regimes with the assumption that the droplets are 
"at rest" with respect to the vapor. As comparison calculations 
show, it is important to use a growth equation designed for the flow 
regime of interest. Otherwise, a serious over-prediction of droplet 
growth may result. The growth equation by Gyarmathy appears to be 
applicable throughout the flow regimes and involves no iteration. 

His expression also avoids the uncertainty associated with selecting a 
mass accommodation coefficient and consequently involves less un- 
certainty in specifying adjustable parameters than many of the other 
growth equations. 


*Now with Martin-Marietta Aerospace,. New Orleans, Louisiana 




INTRODUCTION 


When condensation is encountered in a wind tunnel, the test gas 
may cease to properly simulate flow in ambient-temperature air. 
Understanding the condensation process is essential for maximum 
utilization of cryogenic wind tunnels because the onset of condensa- 
tion effects determines the minimum useful operating temperature which, 
in turn, determines the maximum benefits of cryogenic operation (see 
references 1 to 4 ). An experimental program at the NASA Langley 
Research Center has utilized the 0.3-meter transonic cryogenic tunnel 
to establish a data base for determining the onset of condensation 
effects under a variety of test conditions. Proper analysis of the 
data from the 0.3-m tunnel and extension of the results to larger 
cryogenic tunnels such as the National Transonic Facility (see 
references 5 and 6) require the use of equations for droplet formation 
and growth. With the assumption of no relative motion between the 
droplet and its vapor, this report reviews and evaluates droplet 
growth equations as they apply to flow in transonic cryogenic wind 
tunnels using nitrogen as the test gas. 

The authors would like to express their appreciation of Dr. Gilbert D. 
Stein of Northwestern University for his helpful comments during the 
preparation of this report. 
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heat capacity of liquid, j/(kg.K) 

2 

molecular impingement rate per unit area, l/(m .sec) 

2 

molecular evaporation rate per unit area, 1/ (m .sec) 

Boltzmann's constant, J/K 

Knudsen number, defined in equation (2) 

mean free path, m 

latent heat, J/kg 

molecular mass of condensing vapor 
pressure, N/m^ 

Prandtl number 

distance from the center of a droplet, m 

critical radius of droplet, defined in equation (l), m 

specific gas constant of vapor, j/(kg.K) 

time, sec 

temperature , K 

internal energy, J/kg 

velocity, m/sec 

position, m 

mass accommodation coefficient 
ratio of specific heats 
gas viscosity, kg/ (m. sec) 
thermal conductivity, J/(m.K.sec) 
thermal accommodation coefficient 
density, kg/m 
surface tension, N/m 
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relating to liquid droplet 
relating to saturated value 


relating to vapor 



GROWTH ENVIRONMENT 


Before the droplet growth equations themselves are presented, it is 
useful to review the criterion for droplet growth, the impact of the flow 
regime and the composition of the condensing medium on the growth equations, 
and the assumption of no relative motion between the droplet and its vapor. 


Criterion for Growth 

For a droplet surrounded by its vapor, growth is characterized by the 
vapor molecules condensing on its surface. Continued growth will be ensured 
when the radius of the droplet becomes greater than a critical radius, r^, 
defined in reference 7 as 
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r c ~ R T p,£n(p /p ) 
v v'a v s 


(i) 


where 0 is the surface tension, R^ is the specific gas constant, T^ and p^ 

are the temperature and pressure of the vapor, is the density of the 

droplet, p is the saturation pressure at the vapor temperature, and p /p 
s vs 

is the supersaturation ratio. A droplet with radius r^ = r c will be in 


unstable equilibrium with its vapor while a droplet of radius r^ < r c wi 
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evaporate. The growth of a droplet with r^ > r Q is dependent on the size of 
the droplet relative to the local mean free path and on the composition of the 
surrounding medium. These two factors will be examined individually. 


Flow Regime 


The first factor governing droplet growth is the size of the droplet. 
The Knudsen number, Kn, is a parameter used to define different regimes 
of droplet growth, and is given as 



( 2 ) 
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where £ is the mean free path and 2r^ is the droplet diameter. A growing 
droplet formed by the condensing vapor will initially be smaller than the 
mean free path of the vapor, and therefore Kn » 1. This defines the 
free-molecular regime. The growth of the droplet is determined by molecular 
processes at the surface of the droplet, and kinetic theory provides the 
governing equations in this case. As the droplet grows, it will eventually 
become larger than the mean free path of the surrounding vapor and then can 
be considered to be in a continuum flow regime, where Kn « 1. Macroscopic 
transport phenomena will describe the further growth of the droplet and 
molecular processes will no longer be necessary for the analysis. The 
situation where the droplet diameter is approaching and then surpassing the 
mean free path, Kn — 1, is called the transition regime. Here the importance 
of a molecular description lessens while a macroscopic description becomes 
more important. 

For many condensation problems , it can be assumed that the droplet is 
entirely confined to the free-molecular flow regime during its lifetime, 
thus avoiding the transition and continuum regimes. However, the 
cryogenic nitrogen-gas transonic wind tunnels offer a different environ- 
ment than normally encountered for droplet growth. The two primary differences 
are the higher densities and longer dwell-times involved. In the 0.3-meter 

tunnel for example, where typical stagnation conditions might be T = 120 K 
5 2 

and P = 5 x 10 N/m for a test section Mach number of 0.8, the mean free 

O 

path for nitrogen gas is on the order of 10 m in the test section. For 
many of the past condensation investigations in supersonic and hypersonic 
tests, the mean free path has been on the order of 10 - ^m (see references 8 
and 9). This difference in mean free path implies that the droplets, by 
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equation (2), physically do not have to he as large to attain a given 
value of Kn. Also, since the molecular mean free time between collisions 
is also porportional to mean free path, the smaller mean free path 
will lead to higher collision rates and provide a higher growth rate. 

Furthermore , whereas times for growth in supersonic and hypersonic 
t unn els are often on the order of 10 -1( ' sec (see, for example, reference 8), 

the times for growth in the test section of the 0.3-m transonic tunnel can he as 

_o 

long as 10 sec. When the longer dwe 11-times are coupled with the smaller mean 
free path, the assumption that a free-molecular analysis is sufficient is no 
longer valid. Thus what is needed for transonic cryogenic tunnels is a 
growth equation or equations that describe droplet growth in free-molecular, 
transition, and continuum regimes. 

Composition of Condensing Medium 

While the flow regime is important for droplet growth, the compo- 
sition of the medium surrounding the droplet is also important. The 
0.3-meter tunnel uses nitrogen as the test gas and is consequently a 
single-component system. A more common situation described in the 
literature is the two component system which has the condensing vapor 
dispersed in an inert carrier gas, where "vapor" is used to refer to the 
condensable component and M gas" refers to the noncondensable component. 
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In free-molecular flow, condensation is determined by the kinetics 
of molecular collisions. Once a droplet starts growing, it can only 
continue to grow if the latent heat released by the condensing vapor 
molecules is taken away. Heat can be removed from the droplet in two 
ways - by evaporation or by collisions with either vapor or gas molecules 
where the rebounding molecules receive heat energy from the droplet. The 
presence of an inert carrier gas can increase the probability that a 
vapor molecule will condense on the droplet because the inert molecules 
can be continually removing the latent heat way from the droplet. 

In continuum flow, growth depends on the transport processes of mass 
diffusion and heat conduction. With an inert carrier gas present, heat 
is again carried away by the gas so that vapor- droplet collisions are 
likely to result in condensation. On the other hand, the limiting factor 
for growth with an inert gas present becomes the diffusion of the 
condensable vapor through the gas - the droplet can grow only as fast as 
the vapor molecules can diffuse to it. With no carrier gas present, 
diffusion is not important but the heat generated by condensation must 
be carried away before the droplet can grow. Therefore, with pure 
nitrogen gas, heat conduction is the limiting factor for droplet 
growth in the continuum regime as well as in the free— molecular 
regime . 

Droplet Relative Motion 

An additional assumption used to simplify the problem of droplet 
growth concerns the relative motion of the droplets with respect to the 
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vapor. It is normally assumed that the droplets are "at rest" relative 
to the surrounding medium so that convective heat transfer effects can 
be neglected. Consequently, it is necessary to examine under what 
conditions this "at-rest" assumption will be valid. The breakdown will 
occur when the forces entraining the droplet in the flow are not strong 
enough to keep its velocity the same as the velocity of the vapor. For 
small droplets with Kn » 1, the molecular forces are normally sufficient 
to keep the droplet at the vapor velocity. For a continuum environment 
with Kn « 1, the droplet will stay with the vapor as long as the viscous 
forces tending to keep the velocities the same are as strong as the 
inertial forces on the droplet. The larger the velocity gradient en- 
countered, then the smaller the droplet has to be to stay with the vapor, 
as can be seen by equating the inertial force term to the viscous drag 
term as determined by Stokes equation. Following the reasoning of Wegener 
and Mach in reference 7 s one arrives at the following expression for the 
maximum allowable velocity gradient for negligible velocity difference: 


ax !P d[ ! 



(3) 


where is the velocity of the droplet, w is the velocity of the gas, 
r) is the gas viscosity, is the liquid density of the droplet, and 

r^ is the droplet radius. If it is assumed that the velocity of the droplet 
should not differ from w by more than 1 % , the following table for 
maximum velocity gradients for nitrpgen flow at 100 K can be constructed: 
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r d (m) 

doo/dx (sec”"*") 

10“ 5 

4 x 10° 

10“ 6 

4 x 10 2 

10“ T 

b x 10 1 * 

10- 8 

b x io 6 


Assuming Stokes flow in constructing the above table is a conservative 
assumption because as the Reynolds number of the droplet approaches 
and then surpasses 1, the actual aerodynamic drag on the sphere will 
become larger than that predicted by Stokes flow (see reference 10 ) 
and will consequently allow the droplet to experience larger gradients 
than tabulated above. 

To summarize, for transonic cryogenic wind tunnels using nitrogen as 
the test gas, there is no inert gas present and heat conduction away from 
the droplet will be the limiting factor for the rate of droplet growth. 
Because of the high densities involved and the long dwell-times, droplet 
growth in the continuum regime, as well as the free-molecular regime, must 
be considered. The equations for droplet growth, as presented below, 
will only be valid if the velocity gradients are less than those given 
in the above table for the droplet size of interest. 

EQUATIONS 

Growth equations applicable to the environment of interest in cryogenic 
wind tunnels are outlined along with the fundamental assumptions made by 
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the respective investigators. The equations are presented according to 
their flow regime of applicability. 

Free-Molecular Regime 

Droplets formed by the condensing vapor (homogeneous nucleation) 
normally begin growing in the free-molecular regime, where rarefied gas 
assumptions hold. In particular, kinetic theory (see reference 9 ) pro- 
vides an expression for the mass that crosses unit area per second toward 
one side. This molecular impingement rate per unit area, D, on a droplet 
is given by 


d= — 

(27Tm kT ) ' 
v v 

where the values of the variables are taken from the condensing vapor and 

m is the molecular mass and k is Boltzmann 1 s constant. In order to 
v 

determine the flow rate leaving the droplet, it is assumed that at 
critical radius, r^, the droplet is in quasi^equilibrium and evaporates 
at the same rate as impingement occurs. Then the evaporation rate per 
unit area, D^, can be taken to be equal to the impingement rate on the 
droplet, or 


D 
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K W 

r d a 


(2Trm kTj 1/2 
v d 


( 5 ) 
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where T d is the droplet temperature and P (T ) is the value of pressure 

r d a 

which would just make the droplet of radius r^ a critical-sized droplet. Thus, 
the difference between the condensing, incoming mass and the evaporating, 
outgoing mass is used to calculate the growth rate. In order to complete 
the growth rate equation, a mass accommodation coefficient, a, is also 
needed. The value of a indicates the fraction of molecules that actually 
penetrate into the droplet when they strike it. Altogether, for the 
free-molecular regime the growth rate is given by 


dt 



v 


P (Tj 
d 

27Tm kT , 

v d/ J 


( 6 ) 


A disadvantage with this expression is that the droplet temperature, T^, 
is an unknown quantity. 

reference 11, extends the above results by providing an 
equation to solve for the droplet temperature, T^. He does this by 
balancing the net rate of energy flux to the droplet with the rate of increase 
of its internal energy. The energy equation he obtains, after being 
simplified for a single component gas, is 





+ a(l - f- 
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i T dte: 
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= o 


( 7 ) 
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where 


K 

v 


1 y + 1 

2 Y “ 1 


5 


a is the mass accommodation coefficient, £ is the thermal accommodation 
coefficient, and 


A = (v - Dl 

YRT ’ 

v 

where L is the latent heat of the condensing vapor. The first term in 

equation (j) shows the energy difference between molecules condensing 

on the droplet and molecules evaporating from the surface of the droplet. 

The second term is the energy of reflected vapor molecules from the 

surface of the droplet. The last term is the increase of the internal 

energy of the droplet neglecting the energy needed to raise the 

temperature of the droplet, which is shown to be negligible for the small 

droplets expected in free-molecular flow. In order to solve for r.,, 

d 

T^, and P , an iteration procedure is required. 

In reference 9> Wu has developed a set of equations that bypasses 
the iteration necessary for Hill's expression. Wu starts with kinetic 
theory and defines a new mass accommodation coefficient, a', as 
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and calls it the "effective" mass accommodation coefficient. The last term 
is the ratio of the evaporation rate over the condensation rate. The 
growth equation can then be written as 



Next, Wu specifies the conditions of the droplet and its environment so that 
simplifying assumptions can be made. With these assumptions Wu modifies a' 
by developing four separate equations based on different assumptions. 
Referring to figure 1, the assumed conditions of the droplet include 
thermal and near-thermal equilibrium in addition to saturation and near- 
saturation temperature. For the present case of no inert carrier gas, 

Wu s approximation of the droplet being close to saturation temperature 
is satisfactory. With this assumption, Wu finds that a' becomes 
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vapor pressure of P . (The curve labeled P m figure 1 describes pressures 

v r, 

d 

at which a droplet of radius r^ is just critical-sized. That is, for a given 

value of T^, if P > P^ then the droplet will grow. If P < p , then the 

d r d 

droplet will evaporate . ) 

In another method, Oswatitsch in reference 12 also used kinetic 
theory but in a different way. Oswatitsch took the expression for heat 
transfer and matched it to the expression for mass increase times the 
heat of vaporization. This resulted in the following equation for growth 
rate: 


^d 

dt 






( 11 ) 


This equation was specifically derived for a droplet condensing in its 
vapor without a carrier gas. Oswatitsch also derived an expression 
for droplet temperature by including the Claus ius-Clapeyron vapor-pressure 
curve along with the heat transfer equations. This combination resulted 
in the following equation, 


T - T 

s v 

T - T 

d v 


= i + fi 

2 a 



( 12 ) 
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which may he readily solved for T^. The value of T g used by Oswatitsch 
in equation (12) is just the vapor-pressure curve temperature, T (r = 

S Cl 

P ), and does not take into consideration the droplet radius. Consequently, 
dr d 

T^ and are both independent of droplet size. 

Continuum Regions 

In reference 12, Oswatitsch has also developed an equation for 
continuum flow in which droplet growth is determined by thermal conduction 
on a macroscopic scale. The final expression is 
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2 
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2 A (T, 
v d 


L P; 



(13) 


or 


dr, 

d 

at 


A (T 
v d 


L p. 
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where A^ is the thermal conductivity of the vapor. Unfortunately, Oswatitsch 

does not give an expression for T, for this case as he did for the free- 

a 

molecular case. Consequently, will be assumed in the present paper 
to be equal to T g (r^, ^ v )* Oswatitsch proposed that his first expression. 
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equation (ll), be used while the droplet is in free-molecular flow while the 
second equation (13) > be used for continuum flow, hut there is un- 
certainty in describing droplet growth in the transition regime. 

Fuchs in reference 13 does a survey of evaporation (or condensation 
in reverse) of droplets in the presence of an inert carrier gas. Although 
the inert carrier gas invalidates many of his equations for the present 
single component system, the beginning of Fuchs 1 derivation starts with 
an equivalent to equation (lM before he begins to make assumptions 
appropriate for the presence of an inert carrier gas. 

All Flow Regimes 

Gyarmathy in reference 1^ uses diffusion and heat conduction to 
develop a continuum equation which closely approximates droplet growth in the 
free-molecular regime as well. For a pure vapor Gyarmathy 1 s droplet growth 
equation reduces to 





i + 
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2 y§TT 
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tl.5 Pr 5 Y + 1 
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Kn 



(15) 


where Pr is the Prandtl number for the gas, T is the corresponding 
g s 

value of temperature for a value of on the vapor-pressure curve, and 
Y is the ratio of specific heats. The fact that the applicability of 
this equation extends from the continuum regime to the free-molecular 
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regime makes it the most useful of the equations reviewed. Also, equation 
(15) does not depend on the mass accommodation coefficient. However, 
Gyarmathy's expression neglects (as do Oswatitsch’s continuum equation 
and all of the free-molecular equations) the change in internal energy 
of the droplet due to a varying droplet temperature 

Utt r^ 

% = T- p d c p iT d ’ (16) 


where ‘C is the heat capacity of the liquid. 

While Gyarmathy states that this assumption is excellent for droplets of 

-6 c 

radius 10 m or less, if droplets much larger than 10 -6 m pass through 

significant temperature gradients in the flow, then the droplet growth 

equation should be expanded to include this term. 
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RESULTS AND DISCUSSION 

The droplet growth rates predicted by the various growth rate 

5 2 

equations are presented for a pure nitrogen gas with = 2 x 10 N/m 

and T = 73 K. The different rates are compared to the rate predicted by 

v 

Gyarmathy, equation (15), because equation (15) is designed to be applicable 
in all of the flow regimes. Shaded areas representing uncertainty in the 
equations are also shown to indicate the degree of uncertainty that 
exist in specifying parameters used in the various equations. 

The comparisons between the free-molecular expressions of Hill, 

Wu, and Oswatitsch and the equation of Gyarmathy are shown respectively 
in figures 2, 3» and 4. The uncertainty in what value to assign to 
mass accommodation coefficient, a, is responsible for the large un- 
certainty bands in the three graphs. 

It is not clear what the value of a should be or whether it should 
be the same for one-component or two-component flows. For the 
case of a single component test-gas, steam. Hill found that experimental 
nozzle flow data correlated best with his growth equation if a was 
taken to be 0.04. Similarly, Oswatitsch stated that his own data with 
steam suggested a value for a of 0,03 in order to match the data to 
the predictions of his equations. On the other hand, Wegener and Wu 
in reference 15 recommend a value of a of 0.7 for growth equation developed 
by Wu. The presence of an inert gas in their data may be responsible for 
some of the difference in a. 

In recent years, determination of mass accommodation coefficients 
from film condensation experiments have led to estimations of a much 
higher than 0.0^. For example. Mills and Seban in reference 1 6 have 
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suggested -that the value of ot for steam should be somewhere between 
O.i+5 and unity. This range of values, if also applicable for droplet 
condensation, would seem to support the value recommended by Wegener and 
Wu but would raise some questions about the experimental comparison of 
Oswatitsch and of Hill in addition to questioning the accuracy of 
equation (15 ) for both the free— molecular and transition regimes. 

Consequently, in constructing the uncertainty bands for the free- 
molecular calculations shown in figures 2, 3, and 4, a heavy line is 
drawn for the value of a used by the respective author while the un- 
certainty band extends over a range of a from 0.03 to 0.7. While the 
choices of a vary between Hill, Wu, and Oswatitsch, the uncertainty band 
regions are very similar between the three authors. Also, at i«ii values 
of r, equation (15 ) falls inside the uncertainty regions with the exception 

of the one calculated from Oswatitsch, which does not depend on the radius 

dr 

of the droplet and consequently predicts a constant value of — — . This 

dt 

is in error at the smaller values of r, approaching r , where T = T 

d c d v 

and where the droplet should have a zero growth rate. 

In order to present a comparable uncertainty band in the continuum 
expression of Oswatitsch, an examination of the probable derivation of 
equation (l4) was undertaken to see if any possible difficulty existed 
in specifying the continuum parameters as occurred for a in the free- 
molecular expressions. As shown in the appendix, a problem does appear 
in having to assume that A^ is independent of temperature in order 
to arrive at equation (l4). However, in most cases A^ is nearly 
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porportional to temperature, so if one assumes that 


V T > = w 



( 17 ) 


then the following growth equation can he derived (see the appendix) 


^d 

dt 


VV (T d 


2L p r T 
d d 



(18) 


Nevertheless , as seen in figure 5 where the heavy line represents the 

original equation of Oswatitsch and the uncertainty band extends to the 

line corresponding to equation (l8), there is not a great deal of 

difference between equations (l4) and (l8) for the present P and 

T conditions. Values for the parameters in these equations were taken 

from, or extrapolated from, reference IT. 

As seen in figure 5, equation (l4) and equation (15 ) asymptotically 

approach one another as the continuum regime is reached. This agreement is 

to be expected as the two equations are equal in the limit of Kn + 0 
r 

and — •+ 0. 
r d 

A similar attempt was made to determine what assumptions had 
been made in deriving equation (15). Unfortunately, reference l4 
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cited for Gyarmathy does not contain the necessary information to 
retrace a derivation, so his assumptions could not be examined. In 
particular, it is not clear how the numerical factor in the 

denominator of equation (15) is determined. Nevertheless, as a rough 
approximation to the possible uncertainty in equation ( 15 ), one may 
substitute into the equation a range of values for the parameters 
Pr g and Y (£> the thermal accommodation coefficient, is consistently 
recognized to be unity in the literature so no change was made to this para- 
meter), If the derivation of equation ( 15 ) assumed an ideal gas, then perhaps 
Y should be taken to be l.k instead of the real-gas value of I. 53 . 

Similarly , perhaps Pr should be taken to be the ambient temperature 

o 

value of 0.726 instead of the real-gas value of 0.8U*. The range in 


values that results from the various combinations of the above values 

of y and Pr^ are shown in figure 6, where the heavy line represents 

equation (15 ) with the real-gas values of y and Pr . 

g 

As shown in figures 2 to 6, care must be taken in selecting an 
equation to predict droplet growth rate. Serious over-prediction of 
droplet growth can result if free-molecular expressions are applied 
to either the transistion or continuum regimes or if continuum 
equations are applied to the transistion or free-molecular regimes. 

In the free-molecular equations, setting the value of a, the mass 
accommodation coefficient, is a problem since a range of a between 
0.03 and 0.7 can result in a factor of 3 difference in growth 
rate. Equation ( 15 ) seems to yield reasonable magnitudes of 
growth rate in the free-molecular limit as well as in the continuum 
limit and appears to be well-behaved in the transition regime. Also, 
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equation ( 15 ) has less uncertainty in specifying parameters than the 
free-molecular expressions. Consequently, for the present study equation 
( 15 ) is the preferred representation for droplet growth in all flow 
regimes . 
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SUMMARY OF RESULTS 


It is important to use a droplet growth equation that is designed for 
the flow regime of interest - free molecular, transition, or conti- 
nuum. A serious over-prediction of droplet growth may result if an 
equation is incorrectly applied. The growth equation by Gyarmathy appears 
to be applicable throughout the flow regimes and therefore offers 
an easy solution for droplet growth. The expression by Gyarmathy also avoids 
the uncertainty associated with selecting a mass accommodation coefficient. 
Care, however, must be taken when droplets are larger in radius than 10 ^ m. 
First, possible differences in velocity between the gas and the liquid 
droplet can occur if pressure gradients are encountered because of the 
inertia of the relatively large droplets. Second, the growth rate 
equations have neglected the heat capacity of the droplets and this may 
lead to inaccuracies if the large droplets encounter temperature 
gradients . 
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APPENDIX 


CONTINUUM GROWTH EQUATIONS AND THEIR 
THERMAL CONDUCTIVITY DEPENDENCE 

Expressions for continuum droplet growth may he derived by equating 
the latent heat released during condensation to the heat conduction away 
from the spherical droplet. The exact form of the equation obtained 
depends, however, on whether or not the thermal conductivity, A^, is a 
function of temperature. If A^ is assumed to be a constant, then 
the expression of Oswatitsch, equation (l4), is readily obtained. If A 

v 

is not assumed constant, a modified version of equation (l4) is 
obtained. 

Starting with Fourier* s law of one-dimensional, steady state, heat 
conduction, one may write 


q = -A A 
^ v 


dT 
v 

dr 


(Al) 


where q is the heat flow in J/sec, A is surface area in m^, and 
dT 

v 

is temperature gradient in K/m. For a spherically symmetric problem, 
equation (Al) becomes 


q = - UttX v c 


dT 
\ 

dr 


(A2) 
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Since q is just the total heat flux and is not a function of r or 
T , one may separate variables in equation (A2) and integrate. The 
integration, however, depends upon whether or not A is a function 
of T , and, consequently, r. 

Thermal Conductivity Constant 

If Ay is independent of — and consequently r — the equation 

(A2) may be simply integrated from r = r,tor = °° and from T = T, 

d d 

to T = T . The value of q is simply 

q = UirAy ^ d (T d - Ty) (A3) 


If the heat flow is equated to the latent heat multiplied by the 
rate of mass deposition on the droplet, one may write 



dr 

a 

dt 


p d L 


Combining equations (A3) and (A^) gives 


(AU) 


^d 

dt 


L p d r 


a 

d 



(ih) 
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which is just the expression used by Oswatitsch. 

Thermal Conductivity Variable 

Separating variables and integrating equation (A2) leads to the 
following expression when ^ is a function of temperature 

<A5 ) 


For many gases A (T) can be approximated by 


q = -Inrr. 

^ in 


T 

/ 

V T. 


V T > 


dT 


V T) = x(T a> I: 

d 


(A 6) 


as was done by Penner in reference 18 (see page 282). With the above 
assumption, which appears very adequate for nitrogen, the right side of 
equation (A5) can be integrated to give 


<1 


2 * r a VV 



(AT) 
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Substituting for q. 


dt 



2 L p a r a T a 


(18) 


which is just the corresponding droplet growth expression to equation (ll+) 
assuming is a linear function of temperature. 
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Figure 1.- Conditions of a growing droplet. Point 1 represents thermal 
equilibrium temperature assumption. Point 2 represents 
the saturation temperature assumption for droplet of radius r d . 
(See reference 9.) 



Figure 2.- Free-molecular growth rate of Hill calculated with different values of 
a compared to growth rate predicted by Gyarmathy. 


P = 2 x 10° N/m , T = 73 K, and £ = 1A2 


x 1.0 m. 
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Figure 4.- Free-molecular growth rate of Oswatitsch calculated with different 
values of a compared to growth rate predicted by Gyarmathy. 


P y ■= 2 X 10 5 N/m 2 , T y = 73 K, and l = 1.12 x IQ -8 m. 




r d , radius of droplet, m 


Figure 5.- Continuum growth rate of Oswatitsch assuming thermal conductivity, 
A. v , is either a constant, equation (14), or a linear function of T, 
equation (18), compared to growth rate predicted by Gyarmathy. 

P y = 2 x 10 5 N/m 2 , T y = 73 K, and £ = 1.12 x 10' 8 m. 
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